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Abstract 

Using the Bethe-Salpeter (B-S) equations for heavy baryons Ag, Sq, Hq and 
(Q = b or c), which were estabhshed in previous work, we calculate the 
heavy quark distribution functions in these baryons. The numerical results 
indicate that these distribution functions have an obvious peak at some 
fraction, oq, of the baryon's light-cone "plus" momentum component carried 
by the heavy quark, and that as mq becomes heavier this peak becomes 
sharper and closer to 1. The dependence of the distribution functions on 
various input parameters in the B-S model is also discussed. The results are 
seen to be qualitatively similar to an existing phenomenological model. 



PACS Numbers: ll.lO.St, 14.20.Mr, 14.20.Lq 



I. Introduction 



A significant amount of experimental data liave been accumulated on lepton 
nucleon deep inelastic scattering processes. From these data one can extract in- 
formation on the parton distribution functions in the nucleon, which describe its 
nonperturbative hadronic structure. In comparison with the nucleon, much less 
is known about the parton distribution functions in other baryons such as A, A^, 
and Afe. This is because it is impossible to produce targets of these short lived 
baryons suitable for experiments. Although the parton distribution functions in 
these baryons cannot be studied through deep inelastic scattering processes, it is 
still possible to obtain information on their parton distribution functions by mea- 
suring the fragmentation of quarks to baryons or the decays of these baryons. In 
fact, the parton distribution functions in A, S and A have been calculated in the 
MIT bag model dressed by mesons [||. In this paper, we will study the heavy quark 
distribution functions in heavy baryons. 

The dynamics inside a heavy hadron is simplified by the fact that the light 
degrees of freedom in a heavy hadron are blind to the flavor and spin quantum 
numbers of the heavy quark when its mass is much bigger than the QCD scale 
0. This makes heavy flavor physics a good area for studying nonperturbative 
QCD interactions. In fact, with more measurements on heavy baryons becoming 
available [||, @, ||, , theoretical study of the structure of heavy baryons is becoming 
increasingly important. 

The parton distribution functions are scale-dependent, and their evolution in the 
perturbative region is described by the famous Dokshitzer- Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations [0. However, to determine the parton distribution func- 
tions at some low energy scale (the boundary condition for DGLAP), one needs 
to apply either lattice QCD ^ or nonperturbative effective models. In the heavy 
quark limit, the light degrees of freedom in a heavy baryon have good quantum 
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numbers which can be used to classify heavy baryons. Based on this fact, in a 
previous work we took the heavy baryon to be composed of a heavy quark and a 
hght diquark. With this picture, we estabhshed the B-S equations for the heavy 
baryons Ag and ujq (where uj represents S, S, or and solved these equations 
numerically by assuming that their kernel contains a scalar confinement term and a 
one-gluon-exchange term . It is the purpose of the present paper to calculate the 
heavy quark distribution functions in these heavy baryons in our B-S formalism. 

In Aq and ujq we have 0^ and l"*" diquarks, respectively. Since in the heavy 
quark limit the internal dynamics of a heavy baryon are described by the light 
degrees of freedom, we expect that the heavy quark distribution functions in ujq 
and uJq should be the same. Furthermore, the differences among the heavy quark 
distribution functions in Eg, Sg, and VLq should be caused by SU{?>) flavor breaking 
effects. 

The remainder of this paper is organized as follows. In Section II we derive the 
formulas for the heavy quark distribution functions in various heavy baryons in the 
B-S formalism. In Section III we present numerical results for these distribution 
functions and discuss their dependence on the parameters of the model. We also 
compare our results with the distribution functions of Guo and Kroll [|10]- Finally, 
we give a summary and some suggestions for future work in Section IV. 

II. Formalism for the heavy quark distribution functions 

in Aq, Eq, Hg and VLq 

Based on the picture that a heavy baryon is composed of a heavy quark and a 
light diquark, it was shown that the B-S equation for Ag is 

Xp{v) = Sf{\,P + p) J ^,G{P,p,q)xp{q)So{-X2P + P), (1) 

where Xp{p) is the B-S wave function in momentum space, G{P, p, q) is the kernel, Sp 
and So are the propagators of the heavy quark and light scalar diquark, respectively, 
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with Ai = — — and Ao = — ^ — (Trio and mn are the masses of the heavy quark 
and the hght diquark, respectively), P is the momentum of Aq, and p is the relative 
momentum of the two constituents. 

Similarly, the B-S equation for cuq takes the form 

= Sf(XiP + p) J ^Gp.(P,P,?)Xp(g)'5^'(-A2P + p), (2) 

where Gpj^{P,p, q) is the kernel and 5'^'' is the propagator of the 1+ diquark. 

In the heavy quark limit 

Xop(p) = (/>op(p)maq(^^, s), (3) 

where (f)op{p) is a scalar function, and u\q{v,s) is the Dirac spinor for Aq with 
helicity s and velocity v. 

For cuq there are three scalar functions. A, C and D, in the B-S wave function 

= AB^^iv) + Cv^'pt.B^iv) + Dp'^pt^B^v), (4) 

where B^^v) — :^(7/i+'y/n)75'u('y), andpt = p—{v-p)v. Bj^v) satisfies the constraints 
fB^iy) = B^{y) and v^B^{y) = 0. 

The B-S equations have been solved numerically in the covariant instantaneous 
approximation, assuming the kernels contain a scalar confinement term, Vi, and a 
one-gluon-exchange term, V^, with the following form 

^1 = 77 , - (27r)^(5'(pt - gt) 



3 [(pt-90' + A^'][(Pt-?t)' + Qo] 



^2 - ^771 „ \1 , .^.2ir^„ ° „ \2 , 021' (^) 



where k and a^'^^'' are coupling parameters related to scalar confinement and one- 
gluon-exchange, respectively, where Q\ is a parameter associated with the gluon- 
diquark vertex, and where the parameter \x is introduced to avoid the infra-red 
divergence in numerical calculations, with the limit — > being taken at the end 
of the calculation. 



The twist-2 heavy quark distribution function in = gauge is defined as 



11, 12 



Qia) = V2P+ j '^e-'-''^^- {B\TM^-)l^Mm), (6) 
where ipq is the field operator of the heavy quark Q, = "^(-P'^ + -P^)? 7^ = 
:^(7° + 7^)) i^Q^^') denotes iPq{x) at x"*" = = 0, and \B) represents the heavy 
baryon state with the normahzation {B, P, X\B, P', A') = {2nfPo/mB6x,\'S^{P-P') 
(we have chosen the normahzation convention UbUb = 1). 
We define the two-point function 

M^,(P,fc) = j dx''e-'^^{B\T^Q^{x)^QpmB) (7) 

and then the heavy quark distribution function can be expressed as 

= / ^v^P+5(fc+ - «P+)Tr[7+M(P, k)]. (8) 

The parameter a in Eqs. corresponds to the fraction of the heavy baryon' s 

hght-cone momentum component, P"*", carried by the heavy quark, Q. When it is in 
the range < a < 1, Q{a) measures the probabihty to find the heavy quark with the 
"plus" momentum fraction a. In principle, in a heavy baryon there is a possibility 
to find a heavy antiquark, which is generated from the QCD vacuum. However, 
since we are considering heavy quarks with masses much larger than the QCD scale 
Aqcd, it is very difficult to produce them from the QCD vacuum. Therefore, we 
neglect the heavy antiquark distribution functions. In other words, the valence heavy 
quark distribution function is the same as the heavy quark distribution function, 
Q{a). Since there is only one heavy quark in a heavy baryon we have the following 
normalization condition for Q{a): 

j daQ(a) = 1. (9) 

The two-point function M{P, k) can be evaluated in our B-S framework. After 
some algebra we obtain for Ag in the heavy quark limit the result 

M^^{P,k) = {u^^)^[<Pop{k-X^P)rS^\P~k){u^^)^. (10) 
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The propagator of the hght scalar diquark has the form 



1 

= r— 5— , (11) 

Pi -W^ + ie 



where pi = v ■ p — X2'mAQ and Wp = ypf + mj^. 

The propagator of the heavy quark in the heavy quark hmit has the form 

2{pi + Eo + mD + ie)' ^ ' 

where Eq is the binding energy in the heavy quark hmit. 

We choose to work in the rest frame of Aq, in which we have + fc^ = am/y^ 
(mAg is the mass of Aq) from the 6{k^ — aP^) constraint in Eq. (H). Hence we find 

Sd^ = + 2(1 - a)mA^A;3 + (1 - a)^ml^ - ml]. (13) 

The B-S wave function, (popip), can be expressed in terms of 0op(q'O = / ^0op(p) 
as i: 

hpip) = r— TT- ' w2 u,2^^^ / 7?T3(^i + 2p^^2)0op(gO- (14) 

[pi + Eo + mD + ie){pf - + le) J [2'kY 

The constraint 5{k'^ — aP^) leads to pi = —k^ + (a — l)mAQ, and \pt\'^ = |kp. 

Substituting Eqs. ( p!0i|l4|) into Eq. (^, we integrate out with the aid of 6{k~^ — 

aP^). Furthermore, the component k^ can also be integrated out by choosing the 

appropriate contour. Then we arrive at the following result: 

d^ki 1 



2V2n{l -a) J (27r)2 [Eo + mo + K« - l)^Ag + 2(^-1^ d^^l' + "^d)? 



d^qt 



Viipt - qt) + ({a - l)mAg + i\k^\^ + ml) 

V (a- 1 Wao 



(27r)3 

^2(pt - qt)\ (t>op{qt)\ , (15) 

where \pt\^ = + (kl^^^)^ and where k^^^^ = ^^^^^^^[{a-lfml^-\k^\^ -ml]. 

Substituting Vi and V2 in Eq. (|) into Eq. (0) and integrating out the angular 
coordinates we finally obtain 

Q'Ha) ' rlk,|d|k,| 1 



2v^7r(l ~a)J 2n [Eo + mo + - l)mA^ + ^^^-J^(|kx|2 + ml)]^ 

Q 



~nr 1 — ~, — — r^-T <^op('?i) - (pop[Pt>\ 



3(gg - f^' 
1 



(a - l)mA, 



■(|kj 



+ m 



D) 



X 



2b*l|g*l 



0op(gt) 



where 



m2,]2 + 4(a- l)2|kxPmi^. (17) 



2(1 

In deriving Eq. (|16D we have used the following equations to reduce the three 
dimensional integrations to one dimensional integrations 



and 



Ml? 



4^2 (p2 + ^2 + ^2)2_4p2^2' 



Qtdqt p{qh ^^i\Pt\ + \qt\? + 5'' 



(18) 



(19) 



(27r)3 (p, - g,)2 + 52 J 47r2 2|pt||gi| (|p,| - |g,|)2 + 52 ' 
where p(5'^) is some arbitrary scalar function of qf. 

Now we turn to ujq. The two-point function M{P, k) in this case can be derived 
in a similar way and the result is 

M;^{P, k) = x'aik - \iP)SnlAP - - ^1^)- (20) 

Substituting the B-S equation (0) into Eq. (pO|), using Eqs. (^) and (p!2D, and 
working in the covariant instantaneous approximation, pi = qi (which ensures that 
the B-S equation is still covariant after this approximation), we have 



Tt['j+ M^Q [P, k)] 



1 



Pi + Eq + niD + ie 
d^qt 



[AB" + Cv^pt ■ B + Dptpt ■ 5]7^ 



(27r)^ 



B., 



2pt 



2p? 



+ v^pt ■ B 



Pt 



Pt 



+Pt,Pt ■ B^^P^^^f^^[-CV, + D{V, + 2piV,)] 
^Pt 



(21) 



As for Aq, the B-S wave functions A{pi,pf), C{pi,pf) and D{pi,pf) are related 
to C{pf) and D{pf) through the following equations [^i 

d^qt 



MPi.Pt) 



{pi + Eo + mn + - Wl + it) J (27r)3 

-V2 + Diqt) 



2p? 



2p? 



{Aiq^){V^ + 2piV2) 
iV, + 2piV2)] 



(22) 



m%{pi + Eo + mD + ie){pf - + it) J (27r)3 



{-A{q^)[piV, 



+ {pl + ml)V2]-C{q',) 



2 ^Pt ■ <lt,~. 



Pt 



Vi+pi 



{pt ■ qt)^,~. 



Pi 



+D{q^)^^^[plV, + {p^ + ml)V,] 
pt 



(23) 



DiPi,p') 



ml{pi + Eo + mD + ie){pf - + ie) J (27r)3 
+C{q^) 



{A{q^){V,+piV2) 



Pt-qt , ml{3{pf qt)"^ -pfqt) + '^PtiPf qtYfr 
Pi H TT^ ^2 



Pt 



2pt 



+D{ql) 

, {pt ■ qt) 



mli3{pt ■ Qt)^ - piqf) + 2p'i{pt ■ qtf 



2pf 



{Vi + 2piV2) 



Pt 



-P1V2 



(24) 



and Ai^p'f), C{pt) and D{pl) obey the following three coupled integral equations 

-1 f d^qt 



2Wp{Eo + mn-Wp)J (27r)3 

'"W^ + Diq^) 



{A{q^){V^-2WpV2) 



-C{q^)^' ■ ~ T> i n^„2^ (Pt " qt? ' Pht n~ 



C{p\) 



2pt 



2mlWp{Eo + mn- Wp) J {2tt)^ 
-{{Eo + mD)Wp + ml)V2] 



2p't 



{V, - 2WpV2) 



(25) 



{A{q^)[WpV, 



+C{q^) 



+D{q't) 



Pt 



Pt 



{pt ■ qt? 

Pt 



-WpV, + ((Eo + mn)Wp + m^V^] 



(26) 



Dip', 



2mlWp{Eo + niD - Wp) J (2 



d^qt 



in 



{A{q^){V,-WpV2) 



Pt ■ qt 



WpV,+ 



mli3{pt ■ qt)^ - plq^t) + 2pl{pt ■ qtf ~ 



2pf 



Vo 



-D{ql 
, {Pt ■ Qt) 



Pt ^Ft 
ml{3{pt ■ qtf - pWt) + 2p?(Pt • Qtf 



2pf 



{V, - 2WpV2] 



Pi 



2 WpV2 



(27) 



Once again, we first integrate out with the help of the constraint 6{k~^ — aP~^), 
then we further integrate out by selecting the proper contour which contains the 
pole in k^. With the aid of Eqs. ( ^^ ), and noticing that B^^'^^B^ = — -1= and 
Pt • B'j^pt ■ B = ^^|kp, we obtain: 



|kx|d|kj 



3\/27r(l -a) J 27i Eq + mo - Wp 
+p'tC{p't)f2ip't) +p'Mp't)f3ip't)], 



(28) 



where \pt\ is given in Eq. (|T^, and 

d^qt 



h{p\) 



A{q^,){V, - 2W,y2) - -q'A-C{q^t)V2 + D{ql){V, - 2WpV2 



d^qt 
(27r)3 



Pt Pt 



(29) 
(30) 



d^qt 



A{q^){V,-2WpV2)- 



jPt ■ qtf 

pi 



[-C{q^^)V2 + D{q't){V,-2WpV2) 



(31) 



So far we have been working with the heavy quark limit, mq ^ oo, of the B-S 
equation, but in Ref. we also considered the ^jraq corrections to the B-S equation 
for Aq. To order I/itlq, the heavy quark distribution function in Aq is given by 



Q^'5(a) + AQ^'3( 



a] 



where AQ^'5(a) denotes the ^/mq corrections. To order l/mg, two more scalar 
functions appear in the B-S wave function, which can be related to <f)op{p) in the 
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model of Ref. 0. Consequently, the B-S wave function to order l/mg is given by 

1>t 



Xp{p) = (pop{p) 



1 + 



2mr 



uaJv,s). 



(32) 



With Eq. (|32D we can evaluate the two-point function Mp^{P, k) to order l/mg 
from Eq. (p!0|) . Note that the scalar diquark propagator in Eq. (p^T]) remains 
unchanged to order l/mq. From Eq. after integrating out and k^, we 

have 

Ik^ldlk^l K« - - 2(«-lW (l^^l' + ""^^ 



2V27r(l 



q^dq-, 



a) 



271 [Eo + mB + i(a-l)mAQ + 2(^rTw 

IGlTK 



■(|k^|2 + ml,)]2 



47r^ i + qt + IJ-^r - ^Ptqt 



1 



X 



2bt||gt| 



In 



[a - l)mA + 

^ (a-lW 



-In 





+ 




{\pt 




qt\Y + ^^^ 





+ \qt\? + Ql 


{\Pt 




qt\Y + Ql 



where again \pt\ is given in Eq. (|T^. 

The heavy quark distribution functions have been derived at some scale z/q, which 
is of the order of Aqcd- The QCD running of these functions is controlled by the 
DGLAP equations. In our numerical calculations we will use the evolution code 



provided in Ref. ||T3| to give the heavy quark distribution functions at any higher 
scale. This will be done to the next-to-leading order. Since the QCD interactions 
are flavor independent, we can directly apply this code to the cases of heavy quarks. 
It should be pointed out that the scale z/q cannot be determined in our approach. 
We will treat it as a free parameter and leave it to be determined by the future 
experimental data. 

III. Numerical results 

In this section we will give numerical results for the heavy quark distribution 



(33) 
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functions based on the formulas presented in Section II. The B-S wave functions for 
Aq and ujQ were solved numerically in our previous work by discretizing the inte- 
gration region (0, oo) into n pieces (n is chosen to be sufficiently large and we use 
ra-point Gauss quadrature rule to evaluate the integrals) and solving the eigenvalue 
equation with the kernel Vi and V2 in Eq. (|^) in the covariant instantaneous approxi- 
mation 0. The normalization constants of these B-S wave functions are determined 
by the normalization of Isgur-Wise functions at the zero- recoil point. Substituting 



these numerical solutions into Eqs. (p!6| , p8| , |33D , we can obtain numerical results for 
the heavy quark distribution functions in Aq and ujq. 

In our model we have several parameters, i.e., a^^^\ k, Q^, rriD and Eq. The 
parameter Qq is chosen to be 3.2GeV^ from the data for the electromagnetic form 
factor of the proton |^, The parameters a^'^^^ and k are related to each other 
when we solve the eigenvalue equations with fixed eigenvalues 0. The parameter 
K, is of the order Aqcd'^' where k' is around 0.2GeV^ from the potential model 



15| , p!6| . Thus we let k vary in the region between 0.02GeV and O.lGeV p]. The 
parameters mn and Eq are constrained by the relation rriD + EQ + ^Ei = m\g —rriQ 
for Aq (to order I/itlq) and rriD + Eq = m^^ — rriQ for ujq {rriQ 00). From 
the heavy quark effective theory, mo + Eq and Ei are independent of the fiavor 
of the heavy quark. It was shown that the value of Ei (which is of the order 
Aqcd-^'o) influences the numerical results only slightly. In our numerical calculations 
we use rub = 5.02GeV and rric = 1.58GeV which led to consistent predictions with 
experiments from the B-S equation solutions in the meson case . Consequently we 
have itid + Eq + ^^Ei = 0.62GeV for Aq (where we have neglected 1/ml corrections), 
and mD + Eq = 0.88GeV for Sq, l.OTGeV for Eq, 1.12GeV for VIq in the heavy 
quark limit. The parameter cannot be determined and hence we let it vary 
within some reasonable range. For Aq, we choose niD to be in the range 0.65GeV 
~ 0.75GeV [l^. The axial- vector diquark mass is chosen to vary from O.QGeV to 
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IGeV for Sg, from l.lGeV to 1.2GeV for Sg, and from l.lSGeV to 1.25GeV for fig 
0. With these choices for mo, the binding energy Eq is negative and varies from 
around -30MeV to -130MeV. 



With these parameters, from Eqs. we obtain numerical results for the 

heavy quark distribution functions in Ag and ujq. These results are shown in Figs.l 
to 5. Figs. 1(a), 2(a), 3(a), 4(a) and 5(a) show the dependence on k, for a typical 
rriD, and the dependence on mo for /t=0.06GeV'^. All these results correspond to 
some low energy scale which is of order several hundred MeV. They are evolved to 
some higher scale using the DGLAP equations. In Figs. 1(b), 2(b), 3(b), 4(b) and 
5(b), we show the heavy quark distribution functions at the scale z/^ = lOGeV^, 
where we have taken the low energy scale z/g to be 0.25GeV^, as an example. From 
these figures we can see that for different heavy baryons with the same heavy quark 
flavor the shapes of the heavy quark distribution functions are rather similar. For 
the heavy quark distribution functions at the hadronic scale, there is an obvious 
peak at some "plus" momentum fraction carried by the heavy quark, ao, and this 
peak is much sharper for 6-baryons than c-baryons. Furthermore, ao is much closer 
to 1 for 6-baryons than c-baryons. It can also be seen that QCD evolution makes the 
amplitudes of the peaks much smaller. However, the distinction between 6-quark 
and c-quark distribution functions is still obvious at high u"^. From Figs. 1 to 5 we 
can also see that with the increase of k, which represents the strength of confinement, 
the peaks of the heavy quark distribution functions become lower and their widths 
bigger. This is because the heavy quark behaves more freely when confinement is 
weaker. 



In Ref. Close and Thomas calculated the w-quark distribution function in 
the nucleon in the MIT bag model, where numerical results for au{a) at z/^ = 
5 ~ lOGeV^ were shown for various values of the bag radius. It was shown that 
au{a) has a maximum value around 0.4 ~ 0.5 at a = 0.4 ~ 0.5. In general their 
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argument predicts the peak of the valence distribution (for ground state (L = 0) 
baryons) at {rnH — mD)/mH (where mn is the baryon's mass), which is in the range 
0.87 ~ 0.88 for A;, and 0.67 ~ 0.72 for Ac, respectively. Clearly these expectations 
are in excellent agreement with the peak positions at the hadronic scale, z/q. This 
statement is also true for ujq. We see clearly that with the increase of quark masses, 
the peak positions of the quark distribution functions move closer to 1 and the peak 
values become bigger. 

In Tables 1,2,3,4 for different heavy baryons with different parameters in our 
model, we list the values of ckq, the widths of these peaks, and the total "plus" 
momentum fraction, /, carried by the heavy quark which is defined as 



The width is presented by two values of a at which Q{a) is half of its value at the 
peak, Q{ao). Prom these tables we can see the following points: (i) The values of 
tto, /, and the widths depend not only on the light diquark mass mo but also on 
and mostly the dependence on k is more sensitive for c-baryons. (ii) The value 
of CKo decreases with the increase of m^, which is reasonable, (iii) The widths of 
heavy quark distribution functions for 6-baryons are much smaller than those for 
c-baryons. (iv) For Aq, I/uiq corrections make ckq and / bigger. Hence the heavy 
quark carries more momentum with l/mg corrections included, (v) Numerically, 
in the whole range of parameters we are considering, for A^, ao varies from 0.864 
to 0.890 (0.870 to 0.890) and / varies from 0.800 to 0.857 (0.810 to 0.860) without 
(with) l/rriQ corrections, while for Ac, ao varies from 0.670 to 0.717 (0.690 to 0.760) 
and / varies from 0.522 to 0.651 (0.585 to 0.672); for S^, ao varies from 0.823 to 
0.846 and / varies from 0.809 to 0.834, while for Ec, ao varies from 0.583 to 0.647 
and / varies from 0.493 to 0.589; for S^, ao varies from 0.803 to 0.820 and / varies 
from 0.780 to 0.809, while for Sc, ao varies from 0.503 to 0.580 and / varies from 
0.412 to 0.520; for fib, ao varies from 0.797 to 0.813 and / varies from 0.770 to 0.802, 




(34) 
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Table 1: Values of aQ, widths, and / for A;, with and without l/mg corrections 



mz)=0.65GeV 


K(GeV3) 


0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 




0.890 


0.890 


0.890 


0.890 


width 


0.797 ~ 0.923 


0.833 ~ 0.923 


0.773 ~ 0.940 


0.783 ~ 0.943 


/ 


0.857 


0.860 


0.822 


0.830 


mz5=0.70GeV 


fi:(GeV3) 


0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 




0.876 


0.878 


0.878 


0.884 


width 


0.814 ~ 0.918 


0.816 ~ 0.919 


0.762 ~ 0.935 


0.772 ~ 0.937 


/ 


0.844 


0.848 


0.811 


0.820 


mB=0.75GeV 




0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 


ao 


0.868 


0.870 


0.864 


0.870 


width 


0.789 ~ 0.912 


0.802 ~ 0.914 


0.750 ~ 0.930 


0.762 ~ 0.932 


/ 


0.829 


0.834 


0.800 


0.810 



while for Qc, «o varies from 0.533 to 0.580 and / varies from 0.449 to 0.537. 

As mentioned in Section II, the heavy quark distribution functions should be 
normalized according to Eq. (|^). For some fixed set of parameters in our model, 
the B-S wave functions can be completely determined by solving the B-S equations 
numerically and fixing the normalization of the corresponding Isgur-Wise functions 
at the zero-recoil point. Hence by checking whether the normalization condition 
Eq. (H) is satisfied or not we can test whether our model works well. We have 
performed the integration of Q{a) over a numerically. It is found that for A;, and 
Ac, in the whole range of the parameters, / daQ{a) is very close to 1 (for Ab and Ac 
the deviation from 1 is of order 10~^ and 10^^ ~ 10^^, respectively). Furthermore, 
we find that l/mg corrections make / daQ{a) even closer to 1. Therefore, our 
model for l/mg corrections works in the right direction. For Ub and Uc, although 
/ daQ{a) is not as close to 1 as in the case of Ag, the deviation from 1 still happens 
at order 10~^ at most, which is expected to be improved if l/rriQ corrections are 
also included. Thus the momentum sum rule is satisfied in our model. 
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Table 2: Values of ckq, widths, and / for Ac with and without l/mg corrections 



mD=0.65GeV 




0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 


ao 


0.717 


0.730 


0.710 


0.760 


width 


0.580 ~ 0.807 


0.600 ~ 0.817 


0.443 ~ 0.857 


0.510 ~ 0.870 


/ 


0.651 


0.672 


0.569 


0.625 


mD=0.70GeV 


K(GeV3) 


0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 




0.697 


0.717 


0.697 


0.730 


width 


0.537 ~ 0.797 


0.563 ~ 0.807 


0.417 ~ 0.840 


0.487 ~ 0.856 


/ 


0.617 


0.645 


0.546 


0.605 


mD=0.75GeV 


«;(GeV^) 


0.02 (0(1)) 


0.02 (0(l/mQ)) 


0.10 (0(1)) 


0.10 (0(l/mQ)) 




0.670 


0.690 


0.670 


0.713 


width 


0.500 ~ 0.783 


0.530 ~ 0.797 


0.387 ~ 0.823 


0.550 ~ 0.810 


/ 


0.584 


0.618 


0.522 


0.585 



Table 3: Values of ao, widths, and / for E;, and Ec in the heavy quark limit 



mD=0.90GeV 


K{GeY^) 


0.02 (E,) 


0.10 (Eft) 


0.02 (EJ 


0.10 (Ee) 


ao 


0.846 


0.843 


0.630 


0.647 


width 


0.800 ~ 0.883 


0.750 ~ 0.907 


0.517 ~ 0.720 


0.397 ~ 0.773 


/ 


0.834 


0.834 


0.589 


0.528 


mD=0.95GeV 


K{Gey^) 


0.02 (Eb) 


0.10 (Eft) 


0.02 (E,) 


0.10 (EJ 




0.840 


0.843 


0.610 


0.600 


width 


0.780 ~ 0.880 


0.733 ~ 0.903 


0.470 ~ 0.720 


0.463 ~ 0.723 


/ 


0.820 


0.813 


0.550 


0.493 


mo=1.00GeV 


^•(GeV'^) 


0.02 (E,) 


0.10 (E,,) 


0.02 (E,) 


0.10 (E,) 


ao 


0.830 


0.823 


0.583 


0.597 


width 


0.763 ~ 0.923 


0.740 ~ 0.889 


0.430 ~ 0.707 


0.383 ~ 0.730 


/ 


0.809 


0.814 


0.527 


0.523 
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Table 4: Values of ccq, widths, and / for S;, and Sc in the heavy quark limit 



mzj=1.10GeV 


^(GeV^) 


0.02 (Sb) 


0.10 (Sfc) 


0.02 (5,) 


0.10 (Sc) 


ao 


0.820 


0.817 


0.553 


0.580 


width 


0.773 ~ 0.857 


0.733 ~ 0.877 


0.440 ~ 0.647 


0.343 ~ 0.700 


/ 


0.809 


0.796 


0.520 


0.471 


mD=l.l5GeV 




0.02 (Sb) 


0.10 (Sb) 


0.02 (S,) 


0.10 (SJ 


ao 


0.810 


0.810 


0.533 


0.510 


width 


0.753 ~ 0.850 


0.717 ~ 0.875 


0.393 ~ 0.643 


0.300 ~ 0.690 


/ 


0.795 


0.791 


0.485 


0.440 


mD=1.20GeV 




0.02 (Eb) 


0.10 (5ft) 


0.02 (He) 


0.10 (Se) 




0.803 


0.803 


0.513 


0.503 


width 


0.737 ~ 0.853 


0.700 ~ 0.870 


0.347 - 0.640 


0.257 ~ 0.677 



/ 0.782 0.780 0.445 0.412 



Table 5: Values of ao, widths, and / for and Qc in the heavy quark limit 



mi?=1.15GcV 




0.02 (Qb) 


0.10 (Qb) 


0.03 (Qc) 


0.10 (fic) 




0.813 


0.813 


0.573 


0.580 


width 


0.760 ~ 0.853 


0.727 ~ 0.873 


0.453 ~ 0.667 


0.373 ~ 0.710 


/ 


0.802 


0.772 


0.537 


0. 160 


mD=1.20GeV 


K{Gey^) 


0.02 (Qb) 


0.10 (Qb) 


0.02 (Qe) 


0.10 (fie) 




0.803 


0.810 


0.553 


0.557 


width 


0.750 - 0.894 


0.710 - 0.867 


0.430 - 0.650 


0.343 - 0.697 


/ 


0.786 


0.774 


0.510 


0.464 


mi5=1.25GeV 


^•(GeV'^) 


0.02 m 


0.10 (9.,,) 


0.02 (Q,,) 


0.10 (Q,) 




0.797 


0.797 


0.533 


0.540 


width 


0.733 ~ 0.847 


0.697 ~ 0.863 


0.393 ~ 0.650 


0.313 ~ 0.687 


/ 


0.773 


0.770 


0.474 


0.449 



15 



Now we compare the present results with another phenomenological modeL In 
Ref . ||I0| , also based on the quark-diquark picture, Ag is regarded as composed of a 
heavy quark and a scalar light diquark. The heavy baryon wave function \E'aq {a, k^) 
is proposed as a generalization of the Bauer-Stech-Wirbel [l^ meson wave function 
to the quark-diquark case, 

^^Aq(xi, kx) = NA^a{l - afexp{-b^[kl + ml^{a - af]}, (35) 

where A^Aq is the normalization constant, a = mQ/mp^^, k_|_ is the transverse mo- 
mentum and the parameter h is related to the root of the average square of k_|_. 



(k^). The normalization of the wave function is 

j dad2kx|^AQ(a,kx)p = 1. (36) 
The heavy quark distribution function in this model is defined in the following 

Q{a) = |d2k^|vi>A^(«,k^)|2, (37) 

which leads to 

Q(a) = —^a^{l - aYexp[-2b''ml^{a - af]. (38) 
Here (k^) is defined as 

(kl) = I d«d2k^ki|vI/A^(«,k^)|2, (39) 
from which, using Eq. (^61), we have 

2b^{kl) = 1. (40) 



We choose y (kj_) to vary between 400MeV and 600MeV. For these two pa- 
rameters we plot Q{a) at the hadronic scale in Fig. 6(a) and we plot aQ{a) at 
u"^ = lOGeV^ in Fig. 6(b). We can see from these plots that the heavy quark distri- 
bution functions in this model have a similar shape to those found in the present B-S 
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Table 6: Values of ao, widths, and / for A;, and in the model of Ref. 



v/(ki)(GeV) 


0.4 (A,) 


0.6 (A,) 


0.4 (A,) 


0.6 (A,) 




0.816 


0.764 


0.545 


0.460 


width 


0.766 ~ 0.864 


0.693 ~ 0.830 


0.426 ~ 0.660 


0.306 ~ 0.609 


/ 


0.813 


0.758 


0.540 


0.456 



model. The peak position for Q{a), widths of the distributions, and the total "plus" 
momentum fraction carried by the heavy quark Q, are listed in Table 6. Notice that 
these values are independent of the value of mo in this model. Finally, it can be 
seen from Table 6 that both and / in this model are smaller than those in the 
B-S model. 

IV. Summary and discussion 

Based on our heavy quark and light diquark model for heavy baryons Aq and 
ujQ in the B-S equation approach, we have calculated the heavy quark distribu- 
tion functions in these baryons. This was done in the heavy quark limit for both 
Aq and ujq. Furthermore, l/mq corrections were also included in the case of Aq. 
Our numerical results show that for different heavy baryons with the same heavy 
quark flavor the shapes of the heavy quark distribution functions are similar. At 
the hadronic scale, Aqcd, they have an obvious peak at some "plus" momentum 
fraction carried by the heavy quark, ao, which is much closer to 1 for 6-baryons 
than c-baryons. In addition, the widths of these distributions are much smaller 
for 6-baryons than c-baryons. Furthermore, we have found that the peaks of the 
heavy quark distribution functions become lower and their widths bigger when the 
confinement between the heavy quark and light diquark inside a heavy baryon is 
stronger. The QCD evolution of these distribution functions were also discussed by 
applying DGLAP equations to the next-to-leading order and the results show that 
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the distinction between 6-quark and c-quark distribution functions is still obvious at 
high . The evolution makes the amplitudes of the peaks much smaller, especially 
for the distribution functions in 6-baryons. We also calculated the total "plus" mo- 
mentum fraction carried by the heavy quark, which is much bigger for 6-baryons 
than c-baryons. The dependence of all these numerical results on the parameters in 
our model was discussed in detail. We also checked the normalization condition of 
these distribution functions and found that it is satisfied. In addition, we compared 
our results with an existing phenomenological model for Ag proposed in Ref. |]10| . 



It was found that the shapes of the heavy quark distribution functions in these two 
models are quite similar, although the peak positions in the B-S model are closer to 
1 than those in the model of Ref. |]T0|. We also compared our results for the heavy 
quark distribution functions with the general argument of Close and Thomas based 
on energy and momentum conservation. It was found that the peak positions of the 
quark distribution functions move closer to 1 and the peak values become bigger 
when the quark masses increase. 

Although we cannot determine the exact value of z/q, the hadronic scale at which 
the heavy quark distribution functions are calculated, our results do reflect nonper- 
turbative information in the heavy baryons Aq and loq. We have treated as a 
parameter to be determined by experiment. Even though the heavy quark distribu- 
tion functions cannot be measured directly from deep inelastic scattering processes, 
we can still expect to obtain some information about these functions from the de- 
cays of heavy baryons and fragmentation of heavy quarks to heavy baryons. In 
fact, the latter is directly related to experimental measurements of fragmentation 
processes. However, the theoretical extension to address fragmentation processes is 
more complicated. This is currently under investigation. 



18 



Acknowledgment : 

One of us (Guo) gratefully acknowledges F.-G. Cao for helpful discussions. He 
also thanks O. Leitner and F. Bissey for assistance in preparing the figures. This 
work was supported by the Australian Research Council. 



19 



References 



[1] C. Boros and A.W. Thomas, Phys.Rev. D60 (1999) 074017. 

[2] N. Isgur and M.B. Wise, Phys. Lett. B232 (1989) 113, B237 (1990) 527; H. 
Georgi, Phys. Lett. B264 (1991) 447; see also M. Neubert, Phys. Rep. 245 

(1994) 259 for the review. 

[3] OPAL collaboration, R. Akers et al., Z. Phys. C69 (1996) 195; Phys. Lett. B353 

(1995) 402; OPAL collaboration, K. Ackerstaff et al, Phys. Lett. B426 (1998) 
161. 

[4] UAl Collaboration, C. Albarjar et al, Phys. Lett. B273 (1991) 540. 

[5] CDF Collaboration, F. Abe et al., Phys.Rev. D47 (1993) 2639. 

[6] S.E. Tzmarias, invited talk presented in the 27th International Conference on 
High Energy Physics, Glasgow, July 20-27, 1994; P. Abreu et al, Phys. Lett. 
B374 (1996) 351. 

[7] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; 15 (1972) 
675; G. Altarelh and G. Parisi, Nucl. Phys. B126 (1977) 298; Yu.L. Dokshitzer, 
Zh. Eksp. Teor. Fiz. 73 (1977) 1216 [Sov. Phys. JETP 46 (1977) 641]. 

[8] M. Gockeler et al, Phys. Rev. D53 (1996) 2317; M. Gockeler et al, Nucl. Phys. 
Proc. Suppl. 53 (1997) 81; D. Dolgov et al, Nucl. Phys. Proc. Suppl. 94 (2001) 
303; D. Dolgov, Ph.D. thesis, MIT, Sep. 2000; W. Detmold, W. Melnitchouk, 
J.W. Negele, D.B. Renner and A.W. Thomas, |hep-lat/0103006 . 



[9] X.-H. Guo and T. Muta, Phys. Rev. D54 (1996) 4629, Mod. Phys. Lett. All 
(1996) 1523; X.-H. Guo, A.W. Thomas and A.G. Wilhams, Phys. Rev. D59 
(1999) 116007, Phys. Rev. D61 (2000) 116015. 

[10] X.-H. Guo and P. KroU, Z. Phys. C59 (1993) 567. 

[11] J.C. Collins and D.E. Soper, Nucl. Phys. B194 (1982) 445. 

[12] R.L. Jaffe, Nucl. Phys. B229 (1983) 205. 



20 



[13] M. Miyama and S. Kumano, Comput. Phys. Commun. 94 (1996)185. 

[14] M. Anselmino, P. KroU and B. Pire, Z. Phys. C36 (1987) 89; P. KroU, B. 
Quadder and W. Schweiger, Nucl. Phys. B316 (1988) 373. 

[15] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane and T.-M. Yan, Phys. Rev. 
D17 (1978) 3090. 

[16] H.-Y. Jin, C.-S. Huang and Y.-B. Dai, Z. Phys. C56 (1992) 707; Y.-B. Dai, 
C.-S. Huang and H.-Y. Jin, Z. Phys. C60 (1993) 527; Y.-B. Dai, C.-S. Huang 
and H.-Y. Jin, Z. Phys. C65 (1995) 87. 

[17] F. Close and A.W. Thomas, Phys. Lett. B212 (1988) 227. 

[18] M. Wirbel, B. Stech and M. Bauer, Z. Phys. C29 (1985) 637; M.Bauer, 
B.Stech and M.Wirbel, Z. Phys. C34 (1987) 103. 



21 



Figure Captions 



Fig.l (a) Heavy quark distribution functions at the hadronic scale z/q with the B-S 
equation solved for Ag in the limit mq ^ oo; (b) aQ{a) for Ag at z/^ = lOGeV^ in 
the limit itlq oo. The lines on the right (left) are for A^ (Ac). The solid (dot) lines 
correspond to = O.TOGeV and k = 0.02GeV3 (k = O.lOGeV^). The dashed (dot 
dashed) lines correspond to k = 0.06GeV^ and m/j = 0.65GeV {mo = 0.75GeV). 

Fig. 2 (a) Heavy quark distribution functions at the hadronic scale z/q with the B-S 
equation solved for Ag to order 1/mg; (b) aQ{a) for Ag at z/^ = lOGeV^ to order 
l/rriQ. The lines on the right (left) are for Ab (Ac). The solid (dot) lines correspond 
to ruD = O.TOGeV and k = 0.02GeV=' {k = O.lOGeV^). The dashed (dot dashed) 
lines correspond to k = 0.06GeV'^ and niD = 0.65GeV {mo = 0.75GeV). 

Fig. 3 (a) Heavy quark distribution functions at the hadronic scale z/q with the B-S 
equation solved for Sg in the limit mg oo; (b) aQ{a) for Sg at z/^ = lOGeV^ in 
the limit mg —>■ oo. The lines on the right (left) are for (Sc). The solid (dot) lines 
correspond to mo = 0.95GeV and k = 0.02GeV3 {k = O.lOGeV^). The dashed (dot 
dashed) lines correspond to k = 0.06GeV^ and mo = O.QOGeV {mo = l.OOGeV). 

Fig. 4 (a) Heavy quark distribution functions at the hadronic scale z/q with the B-S 
equation solved for Sg in the hmit mg oo; (b) aQ{a) for Sg at z/^ = lOGeV^ in 
the limit mg oo. The lines on the right (left) are for (He). The solid (dot) lines 
correspond to mo = l.lSGeV and k = 0.02GeV3 (k = O.lOGeV^). The dashed (dot 
dashed) lines correspond to k = 0.06GeV^ and mo = l.lOGeV {mo = 1.20GeV). 

Fig. 5 (a) Heavy quark distribution functions at the hadronic scale z/q with the B-S 
equation solved for Qq in the limit mg oo; (b) aQ{a) for Qq at z/^ = lOGeV^ in 
the limit mg —>■ oo. The lines on the right (left) are for (^c)- The solid (dot) lines 
correspond to mo = 1.20GeV and k = 0.02GeV3 (/t = O.lOGeV^). The dashed (dot 
dashed) lines correspond to k = 0.06GeV^ and mo = l.lSGeV {mo = 1.25GeV). 

Fig. 6 (a) Heavy quark distribution functions for Ag in the model of Ref. []T0| |; 



{h)aQ{a) for Ag at z/^ = lOGeV^ in the model of Ref. The two lines on the 



right (left) are for A;, (Ac). The solid (dashed) lines correspond to y {k±)'^ = 0.4GeV 
{\f{kly = 0.6GeV). 
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Fig.l(b) 
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Fig.2(b) 
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Fig.3(b) 



25 




Fig.4(b) 
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Fig.5(b) 
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